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Abstract New results at four-loop order in perturbative QCD for the first two Taylor coefficients of the 
heavy quark vacuum polarization function are presented. They can be used to perform a precise determi- 
nation of the charm- and bottom-quark mass. Implications for the value of the quark masses are briefiy 
discussed. 
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1 Introduction 

Two-point correlators are of central importance for many 
theoretical and phenomenological investigations in Quan- 
tum Field Theory. As a consequence they are studied in 
great detail in the framework of perturbative calculations. 
Due to simple kinematics (only one external momentum) 
even multi-loop calculations can be performed. The results 
for all physically interesting diagonal and non-diagonal 
correlators including full quark mass dependence are avail- 
able up to O (a^) PEIEJ. 

At four-loop order the two-point correlators can be con- 
sidered in two limits: In the high energy limit massless 
propagators need to be calculated, in the low energy limit 
vacuum integrals ("tadpole diagrams") arise. The evalu- 
ation of the latter in three-loop approximation has been 
pioneered in ref. 0j and automated in ref. 5 . 
Recently first results for physical quantities, which are re- 
lated to four-loop tadpole diagrams, have been obtained. 
The four-loop matching condition for the strong coupling 
constant as at a heavy quark threshold has been calcu- 
lated in ref. [S|[7]. The four- loop QCD contribution to the 
electroweak p-parameter induced by the singlet diagrams 
of the Z-boson self-energy has been computed in ref. jH| ■ 
The detailed knowledge of the heavy quark correlator is 
important for the precise determination of heavy quark 
masses with the help of QCD sum rules. As known from 
ref. P?lll()| . the determination of the charm- and bottom- 
quark mass is further improved, if the four-loop correc- 
tions, hence 0{al), for the lowest Taylor coefficients of 
the vacuum polarization function are available. The sub- 
set of four-loop contributions to the lowest two moments, 
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which involve two internal loops from massive and mass- 
less fermions coupled to gluons, hence of 0{al n'l), has al- 
ready been calculated in ref. jE]. The symbol n/ denotes 
the number of active quark-flavors, contributing through 
fermion-loops inserted into gluon lines. The terms being 
proportional to al nj~ are even known to all orders j in 
perturbative QCD 12 . The symbol ni denotes the num- 
ber of light quarks, considered as massless. In this pa- 
per the complete four-loop contributions originating from 
non-singlet diagrams for the first two Taylor coefficients 
are presented. Singlet contributions have been studied in 
ref. P",TT,TS;. 

In general, the tadpole diagrams encountered in these cal- 
culations contain both massive and massless lines. As is 
well-known, the computation of the four-loop /^-function 
can be reduced to the evaluation of four-loop tadpoles 
composed of completely massive propagators only. Cal- 
culations for this case have been performed in |16II18II17J . 
The outline of this paper is as follows. In section |21 
we briefly introduce the notation and discuss generalities. 
In section Owe discuss the reduction to master integrals, 
describe the solution of the linear system of equations, 
give the result for the lowest two moments and discuss 
briefly the impact on the quark mass determination. Our 
conclusions and a brief summary are given in section 01 



2 Notation and Generalities 



The correlator 11 ^"{q) of two currents is defined as 



W^^iq,]) = I / dxe^''^Q\Tf{x)3%m) , (1) 
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with the current j^{x) — ^{x)^^^^{x) being composed 
out of the heavy quark fields ^{x). The function n^'^{q) 
is conveniently written in the form: 



n^'iq) 
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-q 9^ 



q^q" 



n{q^ 



The vanishing of the longitudinal contribution, as well as 
the confirmation of 11 '^{q^ = 0) = 0, has been used as 
a check of the calculation. The function n{q^) is of phe- 
nomenological interest, because it can be related to the ra- 
tio R{s) = (T(e+ e^ — > hadrons)/cr(e''" e~ — > /i+/i^) with 
the help of dispersion relations: 
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Performing the derivative of eq. lO with respect to q^ one 
obtains on the one hand "experimental" moments 
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ds 
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which can be evaluated from the i?-ratio. On the other 
hand one can define "theoretical" moments 
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which are related to the Taylor coefficients C„ of the vac- 
uum polarization function: 



n{q') 
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(6) 



with z — q^ /{Am ). Symbols carrying a bar denote that 
renormalization has been performed in MS-scheme. The 
Taylor expansion in q^ around q^ = leads to massive tad- 
pole integrals. The first and higher derivatives can be used 
for a precise determination of the charm- and bottom- 
quark mass. However, also the lowest expansion coefficient 
Co has an interesting physical meaning, since it relates 
the coupling of the electromagnetic interaction in differ- 
ent renormalization schemes. 

Sum rules as tool for the determination of the charm- 
and bottom-quark mass have been suggested since long in 
ref. 2H1- This method has then been applied later also to 
the determination of the bottom-quark mass [20] . One of 
the most precise determinations of the charm- and bottom- 
quark mass being based on sum rules in connection with 
the calculation of three-loop moments in perturbative QCD 
has been performed in ref. 9 , with the values of mc(TOc) = 
1.304(27) GeV and rnb{mb) = 4.191(51) GeV as resuhs for 
the charm- and the bottom-quark mass. 

It is convenient to define the expansion of the Taylor 
coefficients C„ of the vacuum polarization function in the 
strong coupling constant as as 



Cr. 



cr+(^)'c!^(^)'cf+(f)^cf+.. 



(7) 



Due to the distinct mass hierarchy of the quarks in the 
Standard-Model, one can consider for the above moments 
one species of quarks as massive {tih — 1) and all lighter 
ones as massless. The number of active quarks n/ is then 
decomposed according to nj = ni + Uh- For convenience 
the symbol nh is kept explicitly in the following. 



3 Calculations and Results 

The Feynman-diagrams have been generated with the help 
of the program QGRAF |2J| • After performing the expansion 
in the external momentum q all integrals can be expressed 
in terms of 55 independent vacuum topologies with addi- 
tional increased powers of propagators and additional ir- 
reducible scalar products. 

In order to reduce this host of integrals to a small set of 
master integrals the traditional Integration-by-parts(IBP) 
method has been used in combination with Laporta's al- 
gorithm 22 , 23 . The resulting system of linear equations 
has been solved with a FORMS |24II25II2()| based program, in 
which partially also ideas described in ref. |23ll27-28, have 
been implemented. The rational functions in the space- 
time dimension d, which arise in this procedure, have been 
simplified with the program FERMAT 2^ . Masking of large 
integral coefficients has been implemented, a strategy also 
adopted in the program AIR [301. In order to achieve the 
reduction to 13 master integrals more than 31 million IBP- 
equations have been generated and solved. This leads to 
integral-tables with solutions for around five million inte- 
grals. Furthermore all symmetries of the 55 independent 
topologies have been taken into account in an automated 
way, by reshuffling the powers of the propagators of a given 
topology in a unique way. Taking into account symmetries 
is important in order to keep the size of the integral-tables 
under control. 

The first four master integrals shown in fig. ^ can be cal- 
culated completely analytically in terms of .T-functions. 
The fifth integral (T52) can be obtained from results of 
ref. mEH. 




Figure 1. Factorized or analytically known master integrals. 
The solid (dashed) lines denote massive (massless) propaga- 
tors. 
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The remaining eight master integrals (fig. [21 have been 
determined with high precision numerics with the differ- 
ence equation method 23 in ref. |32|. Independently they 
have been determined in ref. |33j by constructing an e- 
finite basis. Some of these master integrals have also been 
calculated in ref. PHTT1IHK1IH1 . 




Figure 2. Master integrals where only a few terms of their 
£-expansion are known analytically. The solid (dashed) lines 
denote massive (massless) propagators. 
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Inserting the master integrals and performing renor- 
malization of the strong coupling constant as , the exter- 
nal current and the mass m = to(^) leads to the following 
result (/i — m) : 
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where Riemann's zeta-function C„ and the polylogarithm- 
function Li„(l/2) are defined by 

00 -. 00 -. 

Cn = Er^ a^^d a„ = Li„(l/2) = ^-— . (10) 
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The two numerical constants that appear in eq. and 
eq. were obtained in eq. (19) and eq. (20) of ref. l22j: 



^54,3 — 
262.2 = 



-8445.8046390310298.. 
-4553.4004372195263.. 



and 



(11) 



Numerically the coefficients Co and Ci are given by: 
1.4444 

(1.5863 + 0.1387n,, -I- 0.3714n,) 



Co=p 



-f) (0.0252 n^n^-t- 0.0257 nf- 0.0309 < 
-3.3426 nh - 1.2112 ni + 1.4186) + ... 
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-6.4188 nu - 2.9605 m + 8.2846) + ... . 

(13) 

Using relation jSJ for the first moment at three- and 
four-loop approximations one can assess the influence of 
the new four-loop order on the values of the charm- and 
bottom-quark mass. We first summarize the current status 
of the charm- and bottom-quark masses as obtained in |U] 
from the first moment evaluated to order aj.: 

mdS GeV) = 1.027 ± 0.002 GeV, (14) 

mb(10 GeV) = 3.665 ± 0.005 GeV. (15) 

Note that here and below we display only the uncertain- 
ties coming from the variation of the renormalization scale 
/z in the region /x = 10 ± 5 GeV for the bottom-quark and 
/x = 3 ± 1 GeV for the charm-quark respectively. The ex- 
perimental error is larger and discussed in /9j. 

Let us discuss the influence of the newly computed cor- 
rection on the charm- and bottom-quark masses. In our 
analysis we will closely follow ^Sl. In particular, we will 
borrow from that work the value of the flrst "experimen- 
tal" moment as defined in eq. Q. The latest experimental 
information on R{s) which appeared after publication of 
in] will be taken into account in a future study. 

The inclusion of the four-loop contribution to the func- 
tion Ci leads for the case of the charm-quark to the fol- 
lowing modification of eq. H14|) : 

mc(3 GeV) = 1.023 ± 0.0005 GeV (16) 

For the case of the bottom-quark our result reads: 

TO6(10 GeV) = 3.665 ±0.001 GeV. (17) 

Thus, the four-loop correction does not change the value 
of m{,(10GeV) at all (within our accuracy) but does lead 
to significant decrease of the theoretical uncertainty. 



4 Summary and Conclusion 

The calculation of Taylor expansion coefficients of the vac- 
uum polarization function is important for a precise deter- 
mination of the charm- and bottom-quark mass. Within 
this work we have presented a new result at four-loop or- 
der in perturbative QCD for the first two coefficients of 
the Taylor expansion. For the computation the traditional 
IBP-method in combination with Laporta's algorithm has 
been used in order to reduce all appearing integrals on a 
small set of master integrals. 

With the knowledge of the four-loop contributions the 
theoretical uncertainty is well under control in the view of 
the current and foreseeable precision of experimental data. 
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Note added. 

The results of our calculations as expressed in eq. lO have 
been confirmed in the work |36j . where also a strong re- 
duction of the theoretical error due to the unphysical scale 
dependence has been found. 
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